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Abstract 



Baby Skyrmions are topological solitons in a (2+l)-dimensional field theory which resembles 
the Skyrme model in important respects. We apply some of the techniques and approximations 
commonly used in discussions of the Skyrme model to the dynamics of baby Skyrmions and directly 
test them against numerical simulations. Specifically we study the effect of spin on the shape of 
a single baby Skyrmion, the dependence of the forces between two baby Skyrmions on the baby 
Skyrmions' relative orientation and the forces between two baby Skyrmions when one of them is 
spinning. 



1 Introduction 



The goal of this paper is to study the dynamics of solitons in a (2+l)-dimensional version of 
the Skyrme model. By a soliton we mean a localised, finite-energy solution of a non-linear 
field theory. The Skyrme model is a non- linear field theory for pions in 3+1 dimensions 
with soliton solutions called Skyrmions M. Suitably quantised Skyrmions are models for 
physical baryons. Skyrme's theory is non-integrable and therefore progress in understanding 
Skyrmion dynamics has depended on numerical simulations, approximation schemes or a 
combination of both. This approach has been quite successful in the study of static soliton 
solutions in Skyrme's theory || ||. However, the interactive dynamics of two or more 
Skyrmions, which one needs to understand in order to extract the Skyrme model's predictions 



for the nuclear two-body problem, is more difficult to describe. Certain scattering processes of 
two Skyrmions has been simulated numerically |4j] but the variety of possible initial conditions 
that one could consider is so great that it seems impossible to get an overall picture of the 
scattering behaviour from just a few processes. On the other hand, various approximations 
have been used which typically involve truncating the field theory with infinitely many 
degrees of freedom to a finite-dimensional dynamical system || . Some approximations have 
become widely accepted without, however, having been tested against numerical simulations 
of the full theory. Here we will apply many of the concepts and approximations developed 
in the Skyrme model to our model and directly compare them with numerical simulations. 

Our solitons are exponentially localised in space, a property shared by Skyrmions when 
the physical pion mass is included in the Skyrme model. In our model a soliton has a fixed size 
but arbitrary position and orientation. In two spatial dimensions this corresponds to three 
degrees of freedom, two for the soliton's position and one angle to describe its orientation. 
This should again be compared with Skyrmions, which have a definite size and six degrees 
of freedom, three giving its position in space and three parametrising its orientation. The 
long-range interaction behaviour of our solitons resembles that of Skyrmions in important 
respects: in both cases the asymptotic forces between two solitons depend on their separation 
and their relative orientation and are of the dipole-dipole type. Furthermore, there is a bound 
state of two solitons with a toroidal energy distribution in both cases. When the solitons 
are orientated so that the forces are most attractive and then released from rest they scatter 
through the toroidal configuration and emerge at 90° degrees relative to initial direction of 
their motion. 

Because of all these similarities we call our solitons baby Skyrmions. This term has been 
used quite widely to describe solitons in 2+1 dimensions which resemble Skyrmions in certain 
respects. However, in all the models studied so far the moment of inertia for the rotation 
of a single soliton is infinite so that the rotational degrees of freedom are not dynamically 
relevant. Yet the rotational degrees of freedom are crucial in Skyrmion dynamics. In this 
paper we are therefore particularly interested in those aspects of our model which depend 
on the solitons' orientation. 

After a description of our model and a quick review of its static solutions, discussed 
in detail in |7[], we focus on the following questions. Are there exact solutions of the field 
equations representing spinning baby Skyrmions? How does a baby Skyrmion change its 
shape when it spins? How do the asymptotic forces between two baby Skyrmions depend 
on their relative orientation? What are the forces between two baby Skyrmions when one of 
them is spinning? 

While the questions we investigate here are largely motivated by the (3+l)-dimensional 
Skyrme model we should emphasise that our model is also of interest in the large and growing 
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area of (2+l)-dimensional soliton phenomenology. The solitons in our model have a number 
of properties which are novel in this context. 



2 The Model 



The basic field of our model is a map 

: M 3 h-> Sl o (2.1) 

where M 3 is three-dimensional Minkowski space with the metric diag(l, —1,-1). We set the 
speed of light to 1 and write elements of M 3 as (t, x), where cc is a 2- vector with coordinates 
x\ i=l,2, also sometimes denoted by x and y, We use the notation x a , a = 0, 1,2, to refer 
to both the time and spatial components of (t,x), so t = x°. The target space S^iso is the 
2-sphere of unit radius embedded in euclidean 3-space with the Riemannian metric induced 
by that embedding. Here the suffix 'iso' is used to emphasise the analogy with the target 
space in the Skyrme model, which is often referred to as iso-space. The field (p is a scalar 
field with three components <f) a , a = 1, 2, 3, satisfying the constraint </>-0 = (f)\ + </>| + 0| = 1 
for all x G M 3 . Our Lagrangian density, the static part of which was first considered in || 
and discussed further in R, is 



C = F (^-d a <t>-d a ct) - ^(d Q cj> x drf)- (&*<!> x df>4>) - fi 2 (l - n-<t>)J , 



(2.2) 



where n = (0,0,1) and d a = d/dx a . The constants F, k and fi are free parameters: F 
has the dimension energy and k and /i have the dimension length. It is useful to think of 
F and k as natural units of energy and length respectively and of /i as a second length 
scale in our model. Here we fix our units of energy and length by setting F — k — 1. We 
have already set the speed of light to 1, so we use 'geometric' units in which all physical 
quantities are dimensionless. Thus we cannot set /i to 1 by a choice of units, and we will 
fix its value later after we have discussed its significance. Note also that Planck's constant 
h will be some number, but not necessarily equal to 1. The first term in (|2.2|) is familiar 
from <7-models whose soliton solutions have been studied extensively ||. The second term, 
fourth order in derivatives, is the analogue of the Skyrme term in the usual Skyrme model. 
In we explained in more precise terms in which sense this analogy holds by appealing to 
a general geometric framework for the Skyrme model due to Manton |5). Finally, the last 
term does not contain any derivatives and is often referred to simply as a potential. In three 
(spatial) dimensions the Skyrme term is necessary for the existence of soliton solutions but 
the inclusion of a potential is optional from the mathematical point of view. Physically, 
however, a potential of a certain form is required to give the pions a mass | lOll - By contrast, 
in two dimensions a potential must be included in the above Lagrangian in order to obtain 
soliton solutions. 
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To see this, and to understand the Lagrangian L = J Cd 2 x better, we write it in the 
usual form L = T — V. Here T is the kinetic energy 

T = J ^<b-$+^(<bx xd^cfx, (2.3) 

where the dot denotes differentiation with respect to time, and V is the potential energy 

V = j ^d l cf>-d l (f> + ^(d l c(>xd j( f ) )id j( f ) xd i( f ) )+i2\l-n- ( f>)d 2 x. (2.4) 

We are only interested in fields with finite potential energy and we therefore impose the 
boundary condition 

lim (bit, x) = n (2.5) 
|aj|-*oo 

for all t. As a result we can formally compactify physical space to a 2-sphere S 2 and regard 
the fields <b at a fixed time t as maps from S% to Sf so with an associated integer degree. 
The analytical formula for the degree is 

deg[0] = — / <b-d x <b x d 2 (bd 2 x. (2.6) 

All J 

The degree is a homotopy invariant of the field <b and therefore conserved during time 
evolution. Moreover, it gives a useful lower bound on the potential energy, the Bogomol'nyi 
bound 

V[<b] >47r|deg(0)|. (2.7) 

It is well known that the usual a-model Lagrangian, where both the Skyrme term and 
the potential are omitted, is scale invariant. In particular, this means that a soliton solution 
of the a- mo del can have arbitrary size. The inclusion of the Skyrme term breaks the scale 
invariance and it follows from a simple scaling argument that the potential energy contributed 
by the Skyrme term can be lowered by increasing the scale of a configuration. Thus, in order 
to obtain stable solutions it is necessary to include a potential which, on its own, would 
favour small scales. It follows that soliton solutions in a Lagrangian with both Skyrme term 
and potential have a definite size. 

The precise form of the potential does not matter in such scaling arguments and potentials 



other than the one considered here have been studied in the literature |TTJ . Our potential is 
the analogue of the potential usually chosen in the Skyrme model. It contains a constant \i 
which, in the language of quantum field theory, can be interpreted as the inverse Compton 
wavelength of the mesons in our model. To see this it is best to turn to the equations of 
motion. 

The Euler-Lagrange equations for the Lagrangian L are 

d a (</> x d a cb - d P (b(d P (b-(b x d a <b)) = /i 2 </> x n. (2.8) 
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One simple solution is given by (j>{t, x) = n. It has degree zero and is called the vacuum 
configuration. For a physical interpretation of our Lagrangian it is useful to study the 
equation obeyed by small fluctuation around the vacuum configuration. Decomposing cp 
into a component parallel to the vacuum and a component <p orthogonal to it 

<f> = yjl - p 2 n + p w n + p + 0(ip 2 ), p-n = 0, (2.9) 

one checks that the linearised equation for <p is the massive Klein Gordon equation 

(□+/i 2 )^ = (2.10) 

where □ = <9^<9 M is the wave operator in 2+1 dimensions. In the language of perturbative 
quantum field theory where one quantises small fluctuations around the vacuum, the scalar 
fields pi and p<± therefore correspond to scalar particles or mesons of mass K/j,. We now 
understand \i well enough to fix its value. Our choice is again dictated by the desire to 
reproduce important features of the Skyrme model. There, as in real nuclear physics, the 
Compton wavelength of the pion is of the same order as (in fact, slightly larger than) the 
size of a Skyrmion. Thus, we want to tune \x in our model so that the energy distribution of 
its basic soliton solution is concentrated in a region of diameter « 1///. By trial and error 
we find that this is the case when we set /i 2 = 0.1, which we do for the rest of this paper. 

For later use we note a conservation law that can be read off immediately from the 
equation of motion ( |2.8| ). Taking the scalar product with n on both sides of the equation 
we find that the current 

n</> x d a <t> - (n-^)^f^x d a (p) (2.11) 

has vanishing divergence. The symmetry that leads to this conservation law is 50(2) rota- 
tions of the field <f> about n, which can be written in terms of an angle x as 

(01,02,03) >-> (cosx 0i + sinx0 2 , -sinx0i + cosx02,03)- (2.12) 

We call such a transformation an iso-rotation; the corresponding conserved quantity is called 
isospin and denoted by /: 

I = J n <t> x 4> + (n-di4>)(pi4>-4> x cf))d 2 x. (2.13) 

For a systematic discussion of the symmetries of our model we refer the reader to [[/J, but 
here we need only note that both the Lagrangian L and the degree ( |2.6|) are invariant under 
simultaneous reflections in space and iso-space 

P x ■ (x,y) i-> (-x,y) and (0i, 2 , 3 ) i-> (-0i, 2 , 03)- (2.14) 

Later we will also make use of the invariance of both the degree and the Lagrangian under 
the combination of P x with a with a rotation by 7r in both space and iso-space: 

P y : (x, y) i-> (x, -y) and (0i, 2 , 3 ) i-> (0i, -<f> 2 , fo)- (2.15) 
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3 Static Solutions Revisited 



Time-independent solutions of the equations of motion ( |2.8|) , which are stationary points of 
the potential energy functional V (|2.4j ), were studied in detail in j7|. We briefly recall those 



results which are relevant here. An important class of static solutions of the equations of 
motion consists of fields which are invariant under the group of simultaneous spatial rotations 
by some angle a G [0, 2tt) and iso-rotations by —na, where n is a non-zero integer. Such 
fields are of the form 



<t>{x) 



I sin /(r) cos{n6 — x) 

sin f{r) sin(n# — x) \, (3.1) 
V cos/(r) 



where (r, 9) are polar coordinates in the as-plane and / is function satisfying certain boundary 
conditions to be specified below. The angle x is a l so arbitrary, but fields with different x 
are related by an iso-rotation and therefore degenerate in energy. Thus we concentrate on 
the standard fields where x = 0- Such fields are the analogue of the 'hedgehog' fields in the 
Skyrme model and were also studied in || for different values of \j? . 

The function /, which is called the profile function, has to satisfy 

/(0) = m7i, me Z, (3.2) 

for the field ( |3.1| ) to be regular at the origin; to satisfy the boundary condition ( |2.5| ) we set 

hm/(r)=0. (3.3) 

Here we will only be interested in profile functions where m = 1 (the situation for general 
m is discussed in [J7|). One then finds that the degree of the field fl3.1p is 

deg[0] = n. (3.4) 

For a field of the form ( |3.1|) to be a stationary point of the energy functional V, f has to 
satisfy the Euler-Lagrange equation 

/ n 2 sin 2 A / n 2 sin 2 / n 2 f sin / cos f\ n 2 sin /cos/ 2 
r + / +1 + / r/i sm / = 



(3.5) 



It was shown in that the hedgehog fields (3.1) with n = 1 and n = 2 and profile functions 



satisfying the equation above for those values of n are the absolute minima of V amongst 
all field of degree 1 and 2 respectively. We write (f)^ and <f>^ for those fields in standard 
iso-orientation, i.e. with x = in Q3.1J) , and denote their profile functions by /^ and f( 2 \ 
Translations in physical space and iso-rotations act non-trivially on (p^ and (f)^, so there 
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is a three-dimensional family of minima of the energy functional for both n = 1 and n — 2. 
We call any field obtained by translating and iso-rotating the field 4>^ a baby Skyrmion. 
Baby Skyrmions are the basic solitons of our model; as promised in the introduction they 
have three degrees of freedom: two translational and one rotational. For the total energy, 
or mass, of a baby Skyrmion we find 1.564-47T; the energy density is rotationally symmetric 
and peaked at the baby Skyrmion's centre (where <p^ = (0,0, —1)). The field <p^ (and all 
those obtained by translating and rotating it) may be thought of as a bound state of two 
baby Skyrmions and is described in detail in 0. The energy density is again rotationally 
symmetric but peaked at a distance r ph 1.8 from the centre. The mass is 2.936-47T. 

In subsequent sections we will be interested in the field of a spinning baby Skyrmion and 
the forces between two baby Skyrmions. We recall some simple observation concerning the 
asymptotic behaviour of (jy l > from [7] which are the basis of a remarkably accurate model 
for the dynamical phenomena we will then encounter. For large r, and hence small /, the 
equation ( |3.5| ) for n = 1 simplifies to the modified Bessel equation 

r + -/'-(4+^)/ = °- ( 3 - 6 ) 



A solution of this equation which tends to zero at r = oo is the modified Bessel function 
Ki(fir). Thus, the profile function /W of Q3.5| ) is proportional to K% for large r and we can 
write 

/«(r) ~ P £K x {nr\ (3.7) 

where p is a constant which we will interpret further below. Since the modified Bessel 
function has the asymptotic behaviour 



^W^ e ~T + ° ( ^) (3 - 8) 

the leading term in an asymptotic expansion of f^ 1 ' is proportional to e _/ir / y/r, which shows 
in particular that the (potential) energy distribution of the field <pr* is exponentially lo- 
calised. Most of the analytical results in this paper are based on the observation, made in 
J/J, that the asymptotic field ip^ of <fi^ (defined as in ( |2.9| )) can be interpreted in terms 
of dipole fields. This can be seen as follows. For large r we approximate sin/^ ~ /W and 
cos ~ 1 and, using the asymptotic expression ( ^.7| ) we write the field as 



tpl 1 \x) = ^K 1 (pr)\Bbi{0-x) |- (3-9) 

Or, introducing the orthogonal vectors 

p L = p(cos x, sin x) P 2 = p(-sinx, cosx) (3.10) 



cos(6» - x) 
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and x = x/r, we can write 

<P { a\x) = J^PaXKtifir) = --^vW (H a= 1,2. (3.11) 
However, since the Green function of the static Klein-Gordon equation is K Q (nr), 

(A - /j 2 )Ko(/2r) = -2n5 {2) (x), (3.12) 

we have 

(A - = Pa -V6®(x) a = 1,2. (3.13) 

This equation leads to the interpretation of the asymptotic field as the field produced 
by a doublet of orthogonal dipoles, one for each of the components tpi and , in a linear 
field theory, namely Klein-Gordon theory. The strength of the dipole can be calculated from 
the asymptotic form of One finds, by numerically solving the equation ( 3.5 ), 



p = 24.16. (3.14) 



Once this single number is calculated from the non-linear equation (|3.5|) much can be deduced 



about the dynamics of baby Skyrmions using only the linearised equations of motion. 
4 Spinning Baby Skyrmions 

How does a soliton in two or three dimensions change its shape when it spins? What is the 
interactive dynamics of several solitons when some of them are spinning? From the point of 
view of particle physics these are very natural questions to ask. Yet there are surprisingly 
few soliton models in which these questions has been addressed seriously and even fewer 
in which satisfactory answers have been found. This is partly because the questions do 
not make sense in some of the most popular models. In the much studied abelian Higgs 
model , for example, the soliton solutions, called vortices, are fully characterised by their 
position and have no rotational degrees of freedom. Lumps in the CP 1 model, on the other 
hand, can have an arbitrary orientation, but the moment of inertia associated with changes 
in the orientation is infinite so that the rotational degree of freedom is dynamically frozen 
out. There is a modification of the CP 1 model |13[ in which the solitons, called Q-lumps, 



necessarily spin, but single soliton solutions have infinite energy and in configurations of 
several solitons all solitons have to rotate with the same angular frequency. Thus, the effect 
of relative rotation cannot be investigated. 

In the Skyrme model, spin 1/2 quantum states of a single Skyrmion are models for 
physical nucleons, so the question of spinning Skyrmions has naturally attracted a lot of 
attention. In the first paper on this subject [[14]], it was assumed that a Skyrmion would 
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rotate without changing its shape, and, to obtain the quantum states corresponding to the 
proton, neutron and the A-resonance, it was quantised as a rigid body. Although it was 
quickly pointed out that the classical frequency at which a Skyrmion would have to rotate to 
have spin 1/2 is so large that centrifugal and relativistic effects are important, there seems 
to be no quantitative analysis of these effects in the literature. One reason for this is that, in 
three dimensions, a rotating Skyrmion only has axial symmetry (about the axis of rotation) 
and is not of the 50(3) symmetric hedgehog form of the static solution. Thus, to find the 
exact form of a spinning Skyrmion one needs to solve coupled non-linear partial differential 



equations similar to the ones studied in |T5|| . 

In our two-dimensional model, by contrast, there are solutions describing spinning baby 



Skyrmions which are of the hedgehog form ( |3.1|) . Thus we can study the effect of rotation on 
the soliton's shape, its mass and its moment of inertia by simply solving ordinary differential 
equations. It follows from the 'principle of symmetric criticality' |16| that we can find time- 
dependent solutions of the field equations Q2.8| ) by making the time-dependent hedgehog 
ansatz 



( sin f{r) cos{Q — ut) 
sin f (r) sin(0 - ut) I, (4.1) 



V cos/(r) 

where uj, an arbitrary real number, can be interpreted as the field's angular frequency and 
/ satisfies the boundary conditions ( |3.3| ) and ( |3.2j ) with m = 1. Then the field ( |4.1| ) is 
a solution of the Euler-Lagrange equation (|2.8|) if / satisfies the Euler-Lagrange equation 
obtained from the restriction of the Lagrangian L to fields of the form Q4.ll) . Explicitly this 
is the ordinary differential equation 

(r + (- - u 2 r) sin 2 f)f" + (1 - (uj 2 + 1) sin 2 / + (- - uj 2 r)f sin / cos /)/' 

ry ry £ ry 

— ( uj 2 r) sin / cos / — r/i 2 sin / = 0. (4.2) 

r 

We call solutions of the form ( |4.1| ) spinning baby Skyrmions. The total energy of a spinning 
baby Skyrmion depends on uj and is given by 

M(u) = njr (f 2 + (u 2 + 1) (1 + f' 2 ) sin 2 / + 2/i 2 (l - cos /)) dr. (4.3) 

M(0) is just the mass of a baby Skyrmion calculated earlier and shall henceforth be denoted 
Mq. To study the dependence of M and the field <fi w on u we need to distinguish two regimes, 
uj < fi and uj > \x. This can be seen from the behaviour of the / for large r. In this limit / 
is small and the equation (|4.2|) simplifies to 



f + -f'-(\ + (fi 2 -uJ 2 ))f = 0. (4.4) 



Thus, for uj < fi, f decays exponentially for large r while for uj > \i it is oscillatory. 
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It is interesting to compare the asymptotic field of a spinning baby Skyrmion with the 
field produced by a doublet of scalar dipoles. We then need to consider time-dependent 
dipoles and note here that the equation for scalar fields ip a produced by a rotating pair 
of orthogonal dipoles p a {t), a = 1,2, is the Klein Gordon equation with a time-dependent 
dipole source term: 

(U + ^) Va (t,x) = - Pa (t)-V5^(x). (4.5) 

If the dipoles rotate uniformly at constant angular velocity uj then p a = —uj 2 p a and the 
above equation can be solved explicitly in terms of Bessel functions. We will do this below 
for the two regimes uj < [l and uj > fi and compare the results with the asymptotic form 
of 



4.1 The case uj < fi 



If uj < fi, the equation (^4) is again the modified Bessel equation of first order. As before 



we are interested in a solution / which is exponentially small for large r. Such a solution is 



asymptotically proportional to the modified Bessel function K\{kt), where k = sj — u' 2 , 
i.e. 

f„ZE Kl {Kr) (4.6) 

Z7T 



for some constant p. Then the asymptotic field cp w of the rotating baby Skyrmion (|4.1] j with 
that profile function is 



V f(t 1 x) = ^K 1 ( K r) 




cos^tf - 

sm(9 - ujt) I . (4.7) 



Thus, defining the time-dependent dipole moments 

p ± = p(cos ujt, sin ujt) and p 2 = p(— sinujt, coscut), (4.8) 
the first two components of ^ can be written 

= -^P.(t)-Vffo(«r) a = 1,2. (4.9) 



One checks easily that these fields satisfy the linear equation ( [4.5[ ). Thus, just as in the 
static case, the asymptotic form of the rotating hedgehog solution (|4.1|) may be thought of 
as being produced by a rotating pair of orthogonal dipoles. 

We have solved the radial equation ( |4.2|) for various values of uj < \i. In figure 1 we 
plot the corresponding profile functions. As uj approaches fi from below, the soliton's energy 
distribution becomes more and more spread out, which one may interpret as a centrifugal 
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effect. Note, however, that the initial gradient of the profile functions varies very little and 
that most of the change occurs in the tail, which is well described by the modified Bessel 
function. This is the first indication that one can understand many of the features of a spin- 
ning baby Skyrmion in terms of the field in the asymptotic region where it (approximately) 
obeys linear equations. 

Next we want to understand the dependence of a spinning baby Skyrmion's mass on its 
frequency. Recalling the asymptotic formula ( |3.8|) we see that the energy distribution of a 
baby Skyrmion spinning at u < /i is exponentially localised and that the total mass M(uj) 
is finite. At the critical angular velocity u — fi, however, the equation ( [4.4|) is solved by 
/ = 1/r; the corresponding baby Skyrmion is thus only power-law localised and its mass is 
infinite. For later use we note that the divergent part of the integral in the formula (|4.3| ) is 

7T J r (uj 2 sin 2 / + 2 / u 2 (l - cos /)) dr. (4.10) 

Due to the spreading of the energy density as to approaches /i one needs to integrate to ever 
larger values of r when computing M(uj) numerically using the formula ( ^4.3| ). It is then 
more practical to find constants tq and C such that for r > tq the profile function is well 
approximated by C exp(— kt) / ^pr. Then one integrates the energy density numerically from 
to r and performs the remaining integral analytically, using sin 2 / « / 2 and cos / ~ 1 — f 2 /2. 
We plot the function M(uj) in figure 2. a). It grows rapidly as uj — > \i which is consistent 
with our earlier observation that M(/j) is infinite. 

A further quantity of interest is the conserved charge I discussed at the end of section 
2. For fields of the hedgehog form, where spatial rotations and iso-rotations are equivalent, 
we can interpret I as the angular momentum or spin of the field cf) and we denote it by J in 
this context. One finds that 

J(uj)=uj-2tt J r sin 2 / '(1 + f 2 )dr. (4.11) 

The quantity 

A(w) = J(u)/u (4.12) 

may be interpreted as a moment of inertia. Since a spinning baby Skyrmion changes its shape 
as uj varies, the corresponding moment of inertia changes, too. As mentioned above this effect 
is customarily ignored in the Skyrme model. To check the validity of this approximation in 
our model we define Ao = A(0), for later use. Numerically we find Ao = 2ir- 7.558 

We have calculated J(u) for various values of u < \i using the same technique as described 
above for the computation of M{uS). Like M, J diverges as u — > /x from below, but one checks 
that only the term 

u-2ir [ r sin 2 fdr (4.13) 
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is a divergent integral when lo — //. The divergence is of the same order as that in ( [4.10| ) 
and, comparing coefficients, one is lead to the asymptotic formula 

M-N + fiJ (4.14) 

for some constant N. In figure 2.b) we plot the precise relation between M(u) and J(uo). 
Clearly the graph is well described by the linear formula ( [4.1 4|) already for quite small values 
of u. Note that a linear formula of this form holds exactly for Q-lumps [13|, where the 
constant N is the Q-lumps topological charge. One may interpret it in words as 

"mass of a spinning soliton = constant + meson mass x angular momentum" . 

For small u the mass M{oj) depends quadratically on J{uj) as one would expect for the 
rotation of a rigid body. It is instructive to compare our exact results with the non-relativistic 
rigid body formula 

M = M + (4.15) 

We plot the graph of M(J) is figure 2.b) as well. Note that the rigid body formula is 
only a good approximation to the true mass-spin relation for small spins and small mass 
differences M(u) — M. This observation might be relevant for baryon phenomenology in the 
Skyrme model. There a non-relativistic formula like ( [4.15|) is used to calculate the theoretical 
predictions for the masses of the nucleons and the A particle. However, the nucleon mass 
is about 10% larger than the mass of a Skyrmion and the A is about 40% heavier than a 
Skyrmion. Our calculations indicate that the formula fl4.15| ) is a poor approximation for a 
relative mass difference as large as 40% and that it will generally give too large a value for 
the mass of a spinning soliton at a given angular momentum. 



4.2 The case to > fi 



When uj > n the equation ( |4.4j ) is the (unmodified) Bessel equation of first order, all solutions 
of which are oscillatory for large r. Thus, in terms of k — y/u> 2 — fi 2 we can write the 
asymptotic form of the solution in terms of the Bessel functions of first and second kind 

1 dJ 



Ji(fcr) 



k dr 



(kr) 



nkr 



sin(/cr — -n) 



and 



— — cos(At 7r) 

ixkr 4 



Both these functions may occur, so we write the asymptotic form ip w as 

, ( cos(# — ujt) 



<p u (t, x) = --(pY 1 (kr) + qJ^kr)) 



sin(# — ut) 

V 



(4.16) 



(4.17) 



(4.18) 
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where q and p are constants whose meaning we will explain later. A hedgehog field with the 



asymptotic form (4.18|) has infinite energy (|4.3|) and is thus rather unphysical. Nevertheless 



this solution has a natural interpretation in terms of the dipole model which we will give 
below. 

First, however, we want to study the time evolution of a baby Skyrmion which is given an 
initial angular velocity u > p. To investigate this question we have solved the field equations 
fl2.8|) numerically with initial values (pit = 0) = (fP^ and (pit = 0) = —ujn x <p^\ for uj = 0.5 
and uj = 0.9. The grid for our simulations is a square of 250 x 250 points, extending in both 
the x and y direction from —25 to 25. At the boundary we set the field to the vacuum value 
n and absorb any incident kinetic energy. For both of the initial values of uj we find that 
the baby Skyrmion radiates. In figure 3 we display the field of the baby Skyrmion whose 
initial angular velocity was uj = 0.5, 10 units of time after the start of the simulation. The 
picture clearly shows the sort of spiral pattern which is familiar from dipole radiation in 
linear relativistic field theories [[17j . 

The radiation carries away both energy and angular momentum. As a result the baby 
Skyrmion slows down until the angular velocity has dropped to a value below p in both 
simulations. Figure 4 shows how the total energy for the two simulations decreases with 
time. We have also checked that after 1000 units of time the field has settled down to a 
uniformly rotating field of the form ( f4.1| ) and extract the angular frequency. In the simulation 
where initially uj = 0.5 we now find uj « 0.28 and in the simulation where initially uj = 0.9 
we now find uj ~ 0.3. 

It is instructive to interpret both the solution of the hedgehog form and the numerically 
found solution with the spiral pattern shown in figure 3 in terms of the dipole model. For 
this purpose it is best to combine the asymptotic fields ip 1 and tp 2 into the complex field 
$ = ip l + iip 2 . Then, the Klein-Gordon equation ( |4.5|) , with dipole moments given by ( |4.8| ) 
(where now uj > p), can be written 

(□ + /i 2 )$(t, x) = -pe- iwt {d x + id 2 )S (2 \x). (4.19) 

To solve this equation we separate the time dependence in the form 

$(t,x) = e-™g(x), (4.20) 

so that g has to satisfy the static equation 

(A + k 2 )g(x) = p{d l + td 2 )5 {2 \x), (4.21) 

with k as defined above. Next we need suitable Green functions G of the Helmholtz equation 
in two dimensions, normalised so that 

(A + k 2 )G(kr) = 5 {2 \x). (4.22) 
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A solution which describes an 'outgoing' wave at infinity can be expressed in terms of the 
first Hankel function 



G+(kr) = i^(AT) ~ -i^L-e^-i), (4.23) 

and a solution which describes an 'incoming' wave at infinity is given in terms of the second 
Hankel function 



G-(kr) = l -Hf\kr) ~ i^_L- c -*(^i). (4.24) 
A particular real solution of ( [4.22 ) is thus given by 

\{G + {kr) + G-(kr)) = ^Y (kr), (4.25) 

but to obtain the general solution we should add an arbitrary multiple of the real solution 
of the homogeneous Helmholtz equation 

^{G + (kr) - G-(kr)) = ^J (kr). (4.26) 

Thus, converting to polar coordinates 

and choosing the real Green function ^Y (kr) + ^ J (£r) for some real number q we obtain 
a solution of (|4.19|) 

$ r (t, x) = —(pj^kr) + qY 1 (kr))e i( - e - u ' t \ (4.28) 

which is just the asymptotic form <p u (4.18) of the hedgehog field written in complex notation. 
Thus the hedgehog spinning at uj > \x represents a solution with a radiation field that consists 
of both incoming and outgoing radiation. This is the physical origin of the infinite energy of 
the hedgehog solution. 

It is not difficult to guess which Green function will lead to a solution of (|4.19|) displaying 
the spiral pattern observed in our simulation of the spinning baby Skyrmion. Consider the 
solution constructed from the Green function G + . Using 

Jtr(l) 

!^-(kr) = -kH[ 1) (kr) (4.29) 
dr 

that solution is 



$+(t,aj) = l ^ d -^H^\kr)^p\—e^ kr+e -^ ) . (4.30) 



4 V 8nr 
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Remembering that the real and imaginary part of $ + should be interpreted as the first 
two components of the asymptotic field <p + of a spinning baby Skyrmion we see that, for 
sufficiently large r and a fixed value of t, the direction of cp + is constant along the spirals 
kr = —9. This is precisely the spiral pattern we observed in the numerical simulation of a 
baby Skyrmion spinning with u > fi. The dipole picture shows that it can be accounted for 
in terms of the Green function G + . 

The dipole picture can be used to make sense of many of the qualitative properties of a 
spinning baby Skyrmion. In principle one could also check whether the energy loss through 
radiation plotted in figure 4 can be quantitatively modelled in terms of dipole radiation. How- 
ever, the centrifugal effects in spinning baby Skyrmions change the dipole strength, which 
therefore depends on the baby Skyrmion's angular frequency. This considerably complicates 
the calculations and we therefore have not pursued this path. 

5 The Dipole Model for the Interaction of Baby Skyrmions 

In section 3 we saw that a baby Skyrmion acts like the source of a doublet of scalar dipole 
fields. In [7| it was shown, assuming a certain superposition procedure for well-separated 
baby Skyrmions, that a baby Skyrmion also reacts to the field of a distant baby Skyrmion 
like a doublet of scalar dipoles. There is a similar correspondence between a superposition 
procedure and a linear model for the forces between solitons the Skyrme model. In fl8| 
it was shown that the product ansatz in the Skyrme model (without the pion mass term) 
leads to the same forces between well-separated moving and spinning Skyrmions as a simple 
dipole model for Skyrmions, provided relativistic corrections such as retardation effects are 
included in both approximations. In this section we will ignore such relativistic effects and 
describe the dipole model for slowly moving baby Skyrmions. 

Thus consider two well-separated baby Skyrmions, the first centred at R\ and rotated 
relative to the standard hedgehog (fP^ by an angle xi an d the second centred at R 2 and 
rotated by an angle \i- From |7] we know that, at large separation, the leading term in the 
potential describing the interaction of two baby Skyrmions is the interaction energy of two 
doublets of scalar dipoles in the plane, one situated at R\ and the other at R2 such that 
I i?i — _R 2 | is large compared to The dipole moments of the first dipole are p a , a = 1, 2, 
where 



p 1 =p(cosxi,sinxi) 



p 2 = p(-sinxi,cosxi) 



(5.1) 



and the dipole moments of the second are q a , a = 1, 2, where 



q l = p(cosx2,sinx2) 



<?2 = p(-sfnx2,cosx2). 



(5.2) 
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Then, if the dipoles p 1 and q l and the dipoles p 2 and q 2 interact via a scalar field obeying 
the Klein-Gordon equation with mass /j,, the interaction energy between the doublets p a and 

Q a is 

W = E ^(Pa-V)(g.-V)Ab(Mfl), (5.3) 

a=l,2 

where R = R\ — R2 and i? = Introducing also the relative angle ip = xi — X2 one finds 

2 22 
R) = — cosipAK (fiR) = ULKJfjiR) cos ip. (5.4) 

7T 7T 

In the last step we have used ( |3.12| ) and have omitted the 5^ 2 ^-function term because we are 
only interested in large separations R> 

Thus we obtain the first prediction of the dipole model. The force between two baby 
Skyrmions depends on their relative orientation. In particular two baby Skyrmions in the 
same orientation repel each other; if one is rotated relative to the other by 90° there are no 
static forces; if one is rotated relative to the other by 180° the forces are attractive. The 
forces always act along the line joining the two baby Skyrmions, but in addition there is a 
torque which tends to rotate the relative angle to 180°. In analogy with the terminology 
used in discussing Skyrmion dynamics we call this configuration the most attractive channel. 

To obtain a more quantitative picture we must take into account the mass and the 
moment of inertia of the baby Skyrmions. We assume that the rotations of the individual 
baby Skyrmions are sufficiently slow so that we can approximate the functions M(u) and 
A(o>) by the constants Mq and Ao. Hence our model for the asymptotic dynamics of two 
baby Skyrmions has the Lagrangian 

L dipole = \m R 2 + ^M R 2 2 + h oX l + ^Aoxl ~ R). (5.5) 

In fact the centre of mass position S = (R1+R2) /2 and the angle x — (Xi+X2)/2 are cyclical 
coordinates and decouple from the remaining coordinates. Thus we work in the centre of 
mass frame and set x = and S = 0. Moreover we can introduce polar coordinates (R, <fi) 
for the relative position vector R. Then is also a cyclical coordinate and it is consistent 
to set = = 0. Then we obtain the dynamical system with equations of motion 

1 •■ u? 

-A V> = smipK (/iR) 

l 71 

1 2 3 

-M R = V — cos ipK^fiR). (5.6) 

2 7T 

These equations can easily be solved numerically, and in the next section we will compare 
the soliton trajectories predicted by them with those calculated from the full field equations 
( 12.81) . Some readers may then find it useful to think of the equations (|5.6|) in terms of the 
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coupled motion of a pendulum and a point particle. More precisely the angle ifj may be 
thought of as characterising the angular position of a physical pendulum. Then the first 
equation in ( |5.6| ) specifies the torque acting on the pendulum: it vanishes at the stable 
equilibrium point ip — 7r and the unstable equilibrium point ip = and is maximal when 
ip = 7r/2. Moreover the strength of the torque depends on the 'external' parameter R and 
decreases with increasing R. The second equation in ( [OP can be interpreted as the equation 
for the linear motion of a point particle with position R. The force acting on the particle 
depends on its position, its strength decreasing with increasing R, but is also controlled by 
the 'external' parameter ip. When ip = n (the pendulum's stable equilibrium) the force is 
attractive, tending to decrease R; when if) = (the pendulum's unstable equilibrium) the 
force is repulsive, tending to increase R. 

6 Numerical Simulations 

All the simulations of the field equations to be discussed in this section take place on a 
square grid of 250 x 250 points, extending in both the x and y direction from —25 to 25. The 
initial configurations are constructed from two baby Skyrmion fields using the superposition 
procedure referred to above; thus, these configurations can be characterised by giving the 
individual baby Skyrmions' positions and orientations. The baby Skyrmions' centre of mass 
position and the overall iso-orientation are immaterial for the dynamics, and we take the 
former to be at the grid's origin and usually choose the latter so that the individual baby 
Skyrmions' orientations are equal and opposite. For each simulation we will specify the 
initial relative position and the initial relative orientation, denoted tpQ. We will also consider 
initial conditions where both baby Skyrmions have some non-zero initial velocity. We then 
work in the centre of mass frame, so that the baby Skyrmions' velocities are equal and 
opposite, and we include the effect of Lorentz contraction in our initial configuration. Unless 
specified otherwise we will consider initial velocities along the x-axis, so generically one baby 
Skyrmion is initially in the half plane x > with velocity (—v, 0) and the other in the half 
plane x < with velocity (v, 0), where < v < 1. 

6.1 Scattering from Rest 

Suppose the two baby Skyrmions are initially at rest and well-separated, one centred at (10, 0) 
and the other at (—10, 0). We have calculated the time evolution of the corresponding field 
configuration for a variety of initial orientations ipo G [0,7r] and in figure 5. a) we plot R/2, 
the separation of either baby Skyrmion from the centre, as a function of time. We only show 
the time evolution until the baby Skyrmions collide - the actual collision will be discussed 
in the next section. The qualitative features of the time evolution can easily be understood 
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in terms of the dipole model. 

Consider for example the motion when ipo = it. Then the baby Skyrmions are already 
in the most attractive channel and remain there; hence there is an attractive force between 
their centres throughout their interaction. The plot of the actually observed time evolution 
of R/2 shows exactly such an accelerating motion. When ip Q = n/2 the initial torque is 
maximal, but the initial force vanishes. Thus the relative orientation ip swings through the 
attractive channel ip = n, at which point the force between the baby Skyrmions is maximally 
attractive, but then overshoots and approaches ip = 3tc/2. While ip is close to this value 
the force between the baby Skyrmions is again very small or zero, so that we expect the 
separation parameter R to be a linear function of time here. This is precisely what we see in 
figure 5. a). When ip is decreased further the force between the baby Skyrmions is initially 
repulsive. The baby Skyrmions move apart but at the same time ip increases so that some 
time later the baby Skyrmions are in the most attractive channel. The force is now attractive 
and the baby Skyrmions approach each other again. As the relative orientation oscillates 
the baby Skyrmions experience alternating attractive and repulsive forces and thus perform 
the oscillatory motion most clearly seen in the trajectory for ipo = 0.3 -7T. Finally the baby 
Skyrmions get trapped in the attractive channel and collide. 

When ipQ is decreased further the initial repulsive force increases and as a result the baby 
Skyrmions' separation may initially increase so rapidly that the attractive force which the 
baby Skyrmions experience once they are in the attractive channel is too weak to invert 
their relative velocity. The baby Skyrmions then escape to infinity, which in our simulations 
means that they hit the boundary of the grid. We have investigated boundary effects by 
sending a single baby Skyrmion towards the boundary with velocity v = 0.1 and find weak 
repulsive forces when the baby Skyrmion is approximately 8 units away from the boundary. 
The relevant boundary for the present simulation is at x = ±25, so we interpret simulations 
where R/2 becomes larger than 15 as 'escape to infinity'. This happens for ip < tt/4. In 
our simulations the smallest value of ipo for which the baby Skyrmions ultimately collide is 
ipo = 0.275 -7T, for which the trajectory is also shown in figure 5. a). 

We have also solved the equations ( |5.6| ) numerically for a range of initial values ipo, always 
setting .R(O) = 20, R{0) = ip(0) = 0. Comparing these solutions with the corresponding 
trajectories found in our simulations of the full field equations we find qualitative agreement 
in all cases, but quantitative agreement only for the first part of the trajectories (typically 
< t < 100). Moreover, the dependence of the solutions of ( |5.6|) on ip is exactly as found in 
the field theory. There is a critical value ip c such that for ipo e (ip c , 0], R(t) tends to infinity 
as t — > oo, and for ip e (ip c ,ir], R(t) approaches zero (where the equations are singular) 
after a number of oscillations which becomes arbitrarily large as ipo j ip c . Its numerical value 
is ip c « 0.288545, which should be compared with the value 0.275 found in the field theory. 
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In figure 5.b) we show R and ip as a function of time for ipo = 0.28855 -it. The qualitative 
features of the interplay between the angular motion (ip) and the linear motion (R) discussed 
earlier are clearly illustrated. 



6.2 Head-on Collisions 

The dipole model only describes the asymptotic part of the trajectories discussed so far. We 
have seen, however, that for ip > 0.275 -ir the two baby Skyrmions ultimately adjust their 
relative orientation so that they are in the attractive channel, and collide head on. In the 
next set of simulations we study the outcome of such a head-on collision in the attractive 
channel for a variety of different initial speeds. We fix ip = n and place the baby Skyrmions 
at (7.5,0) and (—7.5,0), giving them initial velocities (— v , 0) and (f,0) respectively, where 
0.1 < v < 0.6. 

In all the simulations the baby Skyrmions merge into the ring-like structure of the 2- 
soliton solution and emerge at right angles to their initial direction of motion. After the 
scattering they move away from each other with their relative orientation still in the at- 
tractive channel. This is the 90° scattering that is now a familiar and apparently generic 
feature of topological soliton dynamics. However, in our model this scattering process is 
accompanied by the emission of a large amount of radiation. In figure 6 we show the energy 
distribution immediately after the collision with v = 0.6: the rings of radiation are clearly 
visible. The radiation carries away so much energy that the baby Skyrmions only escape to 
infinity (for our purposes the boundary of the grid) for v > 0.46. For smaller initial veloci- 
ties the attractive forces between the baby Skyrmions after the collision pull them back and 
they perform another head-on collision, again scattering through 90° and emerging along 
their initial direction of motion and in the attractive channel. The second collision is again 
accompanied by the emission of radiation, so the solitons travel less far than after their first 
collision before they turn round. This process is repeated, but now the motion remains close 
to the ring-like 2-soliton solution at all times. The individual baby Skyrmions are no longer 
distinct, and the motion looks like an oscillatory excitation of the 2-soliton. The emission of 
radiation, however, continues until the kinetic energy has virtually disappeared. The final 
configuration is numerically indistinguishable from the 2-soliton solution. 

The scattering process described above is much more radiative than any observed in 
previous simulations of lump scattering in the CP 1 model M or in other two-dimensional 



versions Skyrme models with potentials different from ours ||11|| . At first sight this is sur- 
prising: the radiation in our model is massive whereas there are massless radiation modes 
in all the comparable models mentioned above. However, our model is also the only one 
in which there are strong attractive forces. Moreover these forces are short-ranged, so that 
the potential energy functional V ( f2~4l ) has a large gradient at configurations consisting of 
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two nearby baby Skyrmions. This means that the time evolution of a configuration in the 
vicinity of those points in the configuration space is not adiabatic. It follows in particular 
that the adiabatic or moduli space approximation proposed for the Skyrme model in || is 
not suitable for describing soliton collisions in our model. 

In the next set of simulations we keep the initial velocity v fixed at 0.5 but vary the initial 
relative orientation ipo between and tt. Initially, the baby Skyrmions are again placed at 
(7.5, 0) and (—7.5, 0). To understand the ensuing interaction processes it is useful to note the 
symmetries of the initial conditions. Recalling that in our conventions the baby Skyrmions 
have equal and opposite initial orientations we first observe that the initial configurations 
are invariant under the reflection P x ( p,14| ). Since P x is a symmetry of the Lagrangian, 
the configurations will remain invariant under that operation during their time evolution. 
Hence, assuming that the baby Skyrmions separate after the collision, we can deduce that 
they must separate along either the x-axis or the y-axis. In other words, if scattering takes 
place, it must be scattering through either 0° (i.e. trivial), 90° or 180°. We also note that, 
if the baby Skyrmions separate along the y-axis, the requirement of invariance under P x 
allows only two possibilities for their individual orientations: either the standard orientation 
where, using the conventions of figure 3, the fields point radially outwards, or the standard 
orientation rotated by tt, where the fields point radially inwards. Thus, after the scattering 
the baby Skyrmions either have the same orientation and are in the most repulsive channel 
or their orientations differ by tt in which case they are in the most attractive channel. 

When ip = or ip = tt the initial configurations are additionally invariant under, 
respectively, the reflection P y (|2.15| ) or the combination of P y with an iso-rotation by tt. 



It follows that the configurations after the interaction must have the same invariances. In 
particular for x/jq = tt we can predict purely on the basis of symmetry that, if the baby 
Skyrmions scatter through 90°, they must be equidistant from the origin after the scattering 
and they must be in the most attractive channel. Of course, this is precisely what we observed 
in our previous simulation. For the other 'special' initial configuration, where tpo = 0, we find, 
however, that the baby Skyrmions scatter through 180°. More precisely they head towards 
each other and slow down until they come to a halt at a separation R « 3. Then they 
turn round and escape to infinity along the line of initial approach. The relative orientation 
ip does not change at all during this process. How, within the constraints imposed by the 
symmetries described above, does the scattering interpolate between 90° scattering and 180° 
scattering as we vary tp from tt to 0? 

In our simulations we find that there is an interval J 90 = [7r/6,7r] such that the baby 
Skyrmions scatter through 90° for ip G J 90 and a smaller interval J 180 = [0, 7r/20] such that 
they scatter through 180 for ip G J 180 . However, as ipo I tt/6 the baby Skyrmions emerge 
from the 90°-scattering with different speeds: the baby Skyrmion moving in the positive 
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y-direction moves faster than the one moving in the negative y-direction. The momentum 
balance is restored by radiation which is emitted predominantly along the negative y -axis. 
This process is presented schematically in figure 7.b). 

When t/j is in the remaining interval J capt = (vr/20, ir/Q) the two baby Skyrmions do not 
separate after the collision but form the oscillatorily excited state of the 2-soliton already 
encountered in the previous set of simulations. However, this time the excited 2-soliton as 
a whole moves in the positive y-direction, and radiation is emitted in the opposite direction 
to restore the momentum balance. This process is sketched in figure 7.c). 

The notion of the most attractive channel is useful for summarising our results. For 
■00 £ I 90 U -f capt the baby Skyrmions merge and scatter through 90°; the closer the two 
baby Skyrmions are initially to being in the most attractive channel the less radiation is 
emitted in the scattering process. For -?/>o £ I 90 the baby Skyrmions therefore escape to 
infinity after the collision; for ipo G J capt they from a coincident configuration which we 
interpret as an asymmetrically deformed 2-soliton solution. The increased radiation and its 
asymmetry is due to that deformation. Finally, for ip e J 180 there is not enough time for the 
baby Skyrmions' relative orientation to adjust before the collision, and as a result the force 
between them is always repulsive. We have not discussed negative values for ipo separately 
here because a scattering process with given negative ipo is related to the processes with — ijj 
by the reflection P y . 

Clearly the precise values of the boundaries of the intervals J 90 , J capt and J 180 depend on 
the baby Skyrmions' initial separation and their speeds. We have repeated the simulations 
just discussed with a different value of the initial velocity v ; the qualitative features discussed 
above and depicted in figure 7 are unaffected by such a change. 

6.3 Scattering with non-zero Impact Parameter 

We have also investigated scattering processes with non-zero impact parameter. We define 
the impact parameter in the usual way as the distance of closest of approach between the 
two soliton centres in the absence of interaction, and denote it by b. 

As in the previous simulations, the nature of the scattering processes to be studied here 
depends crucially on the value of ipo, but, again as before, there are basically two regimes. 
For ifj smaller than some critical value (which depends on the impact parameter and the 
initial velocity) the forces are repulsive and the baby Skyrmions escape to infinity after 
scattering through some non-trivial angle. When ipo is larger than the critical value the 
scattering is more interesting: now the forces are attractive and baby Skyrmions may be 
captured into a bound orbit. The orbiting motion is accompanied by emission of radiation. 
Since the qualitative nature of this process is independent of the precise value of i/jq we have 
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investigated it in detail only for the attractive channel i/jq = n. 

For a quantitative study we have also concentrated on a fixed initial speed, v = 0.4, and 
looked at the dependence of the scattering on the impact parameter b. We already know 
from previous simulations that the the baby Skyrmions will be captured in a bound orbit for 
6 = (the oscillatory excitation of the 2-soliton), but it is clear that the baby Skyrmions will 
escape to infinity after the scattering if b is made sufficiently large. Hence there must be a 
critical impact parameter which separates the two types of scattering. From our simulation 
we conclude that this critical value is approximately 1.5. 

When the baby Skyrmions get captured in a bound orbit their centres trace out ellipse- 
like figures whose perihelion rotates slowly and whose diameter and eccentricity decrease 
as energy is lost through radiation. This is best seen in simulations where the speed v is 
small and the impact parameter is large. In figure 8 we show a part of the baby Skyrmions' 
trajectories for the initial velocity v — 0.1 and the impact parameter b = 12. It is not clear 
from our simulations what the final state of such an orbiting motion is. It is quite possible 
that there is a non-radiating solution of the field equations where two baby Skyrmions orbit 
each other with an angular frequency less than /i. Moreover, our analysis of the spinning 
baby Skyrmion can be extended to 2-solitons, showing that there are finite-energy, spinning 
2-soliton solutions of the hedgehog from provided their angular frequency is less than ji. 
However, contrary to the spinning baby Skyrmion these solutions may well be unstable. In 
our simulations of orbiting baby Skyrmions the kinetic energy decreases all the time but is 
never zero. It is not clear whether either of the periodic but non-radiating solutions described 
above is eventually realised or whether the system will continue to radiate until it reaches 
the static 2-soliton solution. 

6.4 The Effect of Relative Spin 

The last set of simulations we want to discuss addresses the effect of relative spin on the 
interaction of two baby Skyrmions. If one of two well-separated baby Skyrmions is spinning 
and the other at rest the dipole forces between the two average to zero over the spinning 
baby Skyrmion's period of rotation. Thus we expect there to be no net force in this situation. 
This is indeed what we observe: with the stationary baby Skyrmion placed at (—7.5, 0) and 
the other one, spinning at u — 0.2 (with the corresponding profile function) placed at (7.5, 0) 
the relative separation oscillates around the initial value 15 with a small amplitude ~ 0.1 
and angular frequency u. 

Next we investigate the combined effect of relative spin and relative motion of the baby 
Skyrmions' centres. In order to see the influence of spin rather than that of the relative 
orientation on the interaction we should make sure that the interaction processes takes 
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much longer than one period of rotation. Thus we keep the frequency of the spinning baby 
Skyrmion constant at uj = 0.2 and send the two baby Skyrmions towards each other along 
the x-axis with a small speed v. Performing this simulation for v = 0.02 and v = 0.05 and a 
variety of values for we find that the baby Skyrmions always repel each other and scatter 
through 180°. The distance of closest approach is > 9 and the spinning Skyrmion's angular 
frequency does not change during the interaction. This remarkable result can be understood 
in terms of the dipole model as follows. Since the angular velocity is large and the torque 
acting on the relative orientation weak as long as the two baby Skyrmions are not too close 
together we can assume that angular velocity is essentially constant during the interaction 
process. Thus, as the two baby Skyrmions approach each other they experience alternatively 
an attractive and a repulsive force for approximately equal durations given by l/(2a>). This 
leads to a sequence of alternating attractive and repulsive impulses. Since the strength of 
the force between two baby Skyrmions increases with decreasing R, it is clear that every 
attractive impulse will invariably be followed by a stronger repulsive one. However, it may 
happen that at some point a repulsive impulse is strong enough to invert the sign of the 
relative velocity R. Then that impulse will be followed by an attractive impulse which is 
weaker than itself. Thus, for a sufficiently small initial speed v the baby Skyrmions will 
always eventually repel each other and escape to infinity. 

This qualitative explanation in terms of the dipole model can be confirmed by numerical 
solutions of the equations (|5.6| ). In figure 9 we show the trajectories R(t) and ip(t) for a 
particular set of initial values. The plots show clearly that ip indeed remains essentially 
constant during the interaction process, as assumed in our qualitative analysis of the dipole 
model. There is a further interesting effect, though. By carefully measuring R and ip long 
after the interaction we find that, at t > 200, R > v and ip < uj (the difference here is very 
small). Thus rotational kinetic energy has been converted into translational kinetic energy. 

We observe similar effects in our simulations of the full field theory which become more 
pronounced as v is increased. We will describe them for one particular simulation, where 
the baby Skyrmions are again initially at (7.5,0) and (—7.5,0), tp = n/2 and u = 0.2, but 
now v = 0.055. The scattering is again repulsive, but now the distance of closest approach is 
only « 7. At the point of closest approach the spinning baby Skyrmion's angular frequency 
decreases to roughly 0.1 and the other baby Skyrmion begins to spin at that frequency, too. 
The baby Skyrmions then move away from each other, each travelling at the considerably 
increased speed 0.16. One checks that the loss in rotational kinetic energy of ^(0.2) 2 Ao is 
approximately balanced by the gain in translational kinetic energy of (0.16) 2 Mo — (0.055) 2 Mo. 

When v is slightly larger than 0.055 the scattering behaviour becomes sensitive to the 
precise value of the initial relative orientation ip Q . For some values of ipo the baby Skyrmions 
repel as before, but for others their linear momentum is large enough to overcome the 
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repulsive barrier due to their relative spin. Then the baby Skyrmions collide and form the 
oscillating and radiating 2-soliton configuration already encountered in previous simulations, 
but now this configuration rotates as a whole. As in the simulation discussed at the end 
of section 6.3 we are unable to decide on the basis of our numerical results whether all the 
kinetic energy will eventually be radiated away or whether the system will settle down to a 
uniformly rotating 2-soliton solution. 

As v is increased further the baby Skyrmions continue to get captured after the collision, 
even for values of v as large as 0.5. Finally, for v = 0.6, they scatter through « 50° and 
escape to infinity. The collision is again very radiative. Afterwards both baby Skyrmions 
spin at the same angular velocity, whose numerical value is less than 0.1. Thus the radiation 
carries away spin as well as energy. 'Exotic' scattering angles in head-on collisions like 
the one observed here are a familiar feature of lump scattering in the presence of overall 
spin, see || and fl3| . Their occurrence is not surprising. If one or both baby Skyrmions 
are spinning initially, the initial configuration is, in general, no longer invariant under the 
reflection operations discussed in section 6.2. Thus, in the presence of spin, there is no reason 
to expect the scattering angle in a head-on collision to take on special values like 90° or 180°. 



7 Conclusion 



Thinking of baby Skyrmions in terms of pairs of orthogonal dipoles has proven very useful 
for understanding their dynamical properties. The dipole picture reproduces the qualitative 
features of the asymptotic field of a spinning baby Skyrmion, both in the radiative regime 
where the angular frequency uo is larger than /i and in the non-radiative regime u < fi. 
The dipole model is most successful when applied to the interaction between well-separated 
baby Skyrmions. Here it gives an accurate, quantitative description. Moreover it allows one 
to understand in simple qualitative terms the repulsive nature of the forces between baby 
Skyrmions which are spinning rapidly relative to each other and whose centres move slowly 
relative to each other. The dipole model also predicts the possibility of transfer of spin 
and rotational kinetic energy from one baby Skyrmion to the other, in agreement with our 
simulations of the full field theory. This 'spin-exchange' scattering is rather reminiscent of 
the electric charge exchange in dyon scattering |T9"f . 

The dipole picture is based on the linearisation of the non-linear baby Skyrme model. 
When studying the change in shape of a spinning baby Skyrmion or the interaction of 
two baby Skyrmions close together, however, non-linear effects dominate and numerical 
methods become indispensable. Here our simulations also contain a number of lessons for 
the discussion of Skyrmion dynamics. 

We saw that for a rapidly baby spinning Skyrmion centrifugal effects lead to a linear 
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dependence of the mass on the angular momentum. This should be contrasted with the 
quadratic dependence predicted by the non-relativistic, rigid body treatment commonly used 
in the Skyrme model. Our results underline the importance of a more careful treatment of 
centrifugal effects when extracting the baryon spectrum from the Skyrme model. 

The notion of the attractive channel proved useful for understanding a wide range of 
scattering processes. Even though the forces between baby Skyrmions are complicated due 
to their dependence on the relative orientation (as it is for Skyrmions) the torque acting 
on the relative orientation is so strong that interacting baby Skyrmions tend to be in the 
attractive channel by the time they collide. Then they scatter through 90° as if they had 
been in the attractive channel initially. However, this process is the more radiative the 
further the baby Skyrmions were initially from being in the attractive channel. 

Finally, our results illustrate the importance of radiation in soliton dynamics. Our sim- 
ulations show that a particularly large amount of radiation is emitted in a head-on collision 
of two baby Skyrmions. This means that baby Skyrmions only escape to infinity after the 
interaction when the initial relative speed or the impact parameter are quite large. Oth- 
erwise they get captured and, after further emission of radiation, settle down to the static 
2-soliton configuration. We have already pointed out that the radiative nature of this pro- 
cess is probably due to the large gradient of the potential energy functional at configurations 
which describe two baby Skyrmions close together. It may even be possible to relate the 
large amount of radiation accompanying most soliton interactions in our model to the fact 
that the energies of the static soliton solutions are much larger (about 50 %) than the Bo- 
gomol'nyi bound ([2.7]) . In any case, it is clear that the adiabatic approximation to soliton 
dynamics proposed in || would not be useful in our model. Moreover, our results show that 
the inclusion of the pion mass does not necessarily make the adiabatic approximation more 
applicable, as is often claimed. In our model soliton dynamics is much less adiabatic than 
in comparable model with massless mesons. 
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Figure Captions 



Figure 1 

Profile functions for a static baby Skyrmion (bottom curve) and a spinning baby Skyrmion 
with angular frequency (from top to bottom) 0.316,0.3 and 0.2. 



Figure 2 

a) The mass M(u) of a spinning baby Skyrmion in units of 4ir as a function of the angular 
frequency 

b) Mass-spin relationship for a spinning baby Skyrmion. The crosses mark pairs (M, J) 



calculated via ( |4.3| ) and ( |4.11|) for the same value of u; both M and J are plotted in units 
of Air. The solid line is a plot of the function M of J ( 4.15| ); here, too, M and J are plotted 
in units of 4n. 



Figure 3 

Plot of the field of a baby Skyrmion with initial angular frequency uo = 0.5 at time t=10. 
At every lattice site in physical space we plot an arrow of unit length whose direction is 
that of (0i,02) (we identify the 1- and 2-axes in the target space S 2 with those in physical 
space). At the head of the arrow we put a '+' if 03 is positive and a 'x' if 03 is negative. If 
(0i + 0|) < 2 x 10~ 4 no arrow is plotted. Thus the vacuum is represented simply by a '+'. 



Figure 4 

Total energy of spinning baby Skyrmions with initial angular frequencies 0.9 (top) and 0.5 
(bottom) in units of 47r. 



Figure 5 

a) Relative motion of two baby Skyrmions released from rest. R/2 is plotted as a function 
of time for, from top to bottom, ipo = 0.275-7T, 0.3-7T, 0.4-7T, 0.5-7T and -00 = 7t "■ 

b) Prediction of the dipole model for the relative motion of two baby Skyrmions released 
from rest. R/2 and ip as a function of time for i/jq = 0.28855-7T. 
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Figure 6 

Total energy density shortly after the head-on collision of two baby Skyrmions with initial 
speeds v = 0.6. 

Figure 7 

Head-on Collisions 

a) Sketch of baby Skyrmion velocities and radiation emitted shortly after scattering with 

^0 = 7T 

b) Sketch of baby Skyrmion velocities and radiation emitted shortly after scattering with 
^o = | 

c) Sketch of 2-soliton velocity and radiation emitted shortly after scattering with -?/>o £ I capt 

d) Sketch of baby Skyrmion velocities shortly after scattering with ^0 = 

Figure 8 

Trajectories of two baby Skyrmions with initial speed v — 0.1 and impact parameter b = 12 

Figure 9 

The effect of spin: trajectories R(t)/2 and ip(t)/2n calculated from the dipole model with 
initial values R(0) = 15, R(0) = 0.06, ip = 0.5-vr and ^(0) = oo = 0.2 
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